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In 2012, a novel betacoronavirus, designated Middle East re- 
spiratory syndrome coronavirus or MERS-CoV and associated with 
severe respiratory disease in humans, emerged in the Arabian 
Peninsula. To date, 108 human cases have been reported, including 
cases of human-to-human transmission. The availability of an 
animal disease model is essential for understanding pathogenesis 
and developing effective countermeasures. Upon a combination of 
intratracheal, ocular, oral, and intranasal inoculation with 7 x 10° 
50% tissue culture infectious dose of the MERS-CoV isolate HCoV- 
EMC/2012, rhesus macaques developed a transient lower respira- 
tory tract infection. Clinical signs, virus shedding, virus replication 
in respiratory tissues, gene expression, and cytokine and chemo- 
kine profiles peaked early in infection and decreased over time. 
MERS-CoV caused a multifocal, mild to marked interstitial pneu- 
monia, with virus replication occurring mainly in alveolar pneumo- 
cytes. This tropism of MERS-CoV for the lower respiratory tract 
may explain the severity of the disease observed in humans and 
the, up to now, limited human-to-human transmission. 
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[I June of 2012, a novel betacoronavirus, associated with severe 
respiratory disease in humans emerged in the Middle East (1, 
2), which is closely related to betacoronaviruses circulating in 
bats (3, 4). The first isolate of Middle East respiratory corona- 
virus (MERS-CoV) (5), HCoV-EMC/2012, was obtained from 
a patient with a fatal pneumonia and acute renal failure. To date, 
107 additional human cases have been identified, of which 49 
were fatal (6). Aside from cases in Saudi Arabia, Qatar, Jordan, 
and the United Arab Emirates, imported cases have been iden- 
tified in the United Kingdom, Germany, France, Tunisia, and 
Italy (6). Although no information is available on the source or 
route of primary transmission of MERS-CoV, human-to-human 
transmission has been recorded (7-9). Clinical data on human 
cases of MERS-CoV infection are currently sparse, but it 
appears that this virus mainly causes severe lower respiratory 
tract disease, occasionally accompanied by renal disease. The 
severity of disease distinguishes MERS-CoV from other coro- 
naviruses circulating in the human population, HCoV-229E, 
HCoV-0C43, HCoV-NL63, and HCoV-HKU1, which are gen- 
erally associated with upper respiratory tract infections. Instead, 
MERS-CoV appears to be more similar to the severe respiratory 
disease caused by severe acute respiratory syndrome (SARS)- 
CoV. 

In vitro studies have shown that MERS-CoV replicates effi- 
ciently in nonciliated cells in the primary human airway epithe- 
lium (10), and in ex vivo human lung cultures MERS-CoV 
replicated in bronchial, bronchiolar, and alveolar epithelial cells 
(11), in line with the observed respiratory disease in humans. The 
recently defined receptor for MERS-CoV, dipeptidylpeptidase 4 
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(DPP4), is generally expressed in endothelial and epithelial cells 
and has been shown to be present on cultured human nonciliated 
bronchiolar epithelium cells (12), providing further information 
on the respiratory tropism of MERS-CoV. 

Animal models that recapitulate human disease are essential 
for understanding pathologic processes involved in disease pro- 
gression. Moreover, these models are instrumental for the de- 
velopment of prophylactic and therapeutic countermeasures. We 
have previously shown that rhesus macaques inoculated with 
a high dose of MERS-CoV isolate HCoV-EMC/2012 developed 
a respiratory disease reminiscent of that observed in humans 
(13). To increase our understanding of the pathogenesis of 
MERS-CoV in the absence of clinical and pathological data 
from human patients, we present herein a more detailed analysis 
of the extent of virus replication, the histopathological changes 
in the respiratory tract and changes in systemic (peripheral blood 
mononuclear cell, PBMC) and local (lung tissue) gene expres- 
sion of MERS-CoV-infected rhesus macaques. 


Results 


Clinical Signs and Hematological Changes in Animals Inoculated with 
MERS-CoV. Six rhesus macaques were inoculated with 7 x 10° 


Significance 


The Middle East respiratory syndrome coronavirus (MERS-CoV) 
is the latest emerged coronavirus causing severe respiratory 
disease with a high case fatality rate in humans. To better 
understand the disease caused by MERS-CoV, we developed 
a rhesus macaque model. Infection of rhesus macaques with 
MERS-CoV resulted in the rapid development of a transient 
pneumonia, with MERS-CoV replication largely restricted to the 
lower respiratory tract. This affinity of MERS-CoV for the lungs 
partly explains the severity of the disease observed in humans. 
The MERS-CoV rhesus macaque model will be instrumental in 
developing and testing vaccine and treatment options for an 
emerging viral pathogen with pandemic potential. 
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Hematological changes in rhesus macaques inoculated with MERS-CoV. Total number of white blood cells (Left), neutrophils (Center), and lym- 


phocytes (Right) were determined in blood samples obtained from animals 0, 1, 3, and 6 dpi. Of note, CoV-2, CoV-4, and CoV-6 were euthanized 3 dpi, and 
thus no samples are available for 6 dpi. Each line represents one animal: CoV-1 (@), CoV-2 (), CoV-3 (A), CoV-4 (W), CoV-5 (4), CoV-6 (C). 


TCIDs 9 (50% tissue culture infectious dose) of MERS-CoV 
isolate HCoV-EMC/2012 (13). Clinical signs upon inoculation 
with HCoV-EMC/2012 were transient and varied between ani- 
mals, as specified in Table S1. Analysis of blood samples col- 
lected 0, 1, 3, and 6 d postinfection (dpi) indicated increased 
white blood cell counts in five of six animals on 1 dpi (Fig. 1). 
Further analysis indicated that this was the result of an increase 
in the number of neutrophils. Moreover, there was a reduction in 
the number of lymphocytes. By 3 dpi, blood cell counts had returned 
to numbers comparable to those before inoculation (Fig. 1). 


Viral RNA in Swabs and Bronchoalveolar Lavages. Nasal, oropha- 
ryngeal, urogenital, and rectal swabs were obtained from the 
animals 1, 3, and 6 dpi. On 1 and 3 dpi, nasal swabs from all 
animals were positive for the presence of viral RNA by quanti- 
tative RT-PCR (qRT-PCR); only one animal had a positive nasal 
swab on day 6 postinfection (Fig. 24). Oropharyngeal swabs 
obtained from CoV-3 on days 3 and 6 postinfection and obtained 
from CoV-6 on day 6 postinfection were positive. All urogenital 
and rectal swabs were negative with the exception of a urogenital 
swab collected from CoV-2 on day 1 postinfection. Bronchoalveolar 
lavages (BAL) were also performed 1, 3, and 6 dpi and analyzed 
for the presence of viral RNA. BAL obtained 1 dpi were pos- 
itive by (RT-PCR in five of six animals; 3 dpi three of six BALs 
were positive, and 6 dpi two of three BALs contained viral 
RNA (Fig. 2B). 


Widespread Dissemination of Virus Through the Respiratory Tract. 
We have previously reported the presence of viral RNA and 
infectious virus throughout the lungs of inoculated macaques, 
with viral load decreasing between 3 and 6 dpi (13). In addition 
to the different lung lobes, we have previously analyzed several 
other tissues, including tissues of the upper respiratory tract, 
lung lesions, and kidney (for complete list, see Materials and 
Methods). Although there was some variation between the 
different macaques with regard to the presence of viral RNA, 
we could consistently detect the virus by FRT-PCR in the nasal 
mucosa, trachea, and mediastinal lymph nodes on 3 dpi (Fig. 3). 
Furthermore, we could detect viral RNA in conjunctiva, tonsils, 
oronasopharynx, and in the left and right bronchus. Viral loads 
were lower in these tissues by 6 dpi and viral RNA could no 
longer be detected in the nasal mucosa or conjunctiva at that 
time (Fig. 3). Viral RNA could not be detected in kidney or 
bladder tissue samples. 

To further prove that virus replication occurred in the lungs of 
MERS-CoV-inoculated rhesus macaques, transmission electron 
microscopy was performed on lung lesion samples collected 6 
dpi. At this time point, coronavirus particles could be detected in 
alveolar pneumocytes (Fig. 4). Moreover, qRT-PCR analysis of 
RNA extracted from the right lower lung lobe and lung lesions 
obtained 3 and 6 dpi, indicated the transcription of subgenomic 
mRNA (ORFS) as well as ORF1 RNA in samples collected from 
all animals and at both time points (Fig. S1), another indicator of 
active virus replication in the lower respiratory tract. 
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Pathology of MERS-CoV in the Respiratory Tract. Upon necropsy of 
animals 3 and 6 dpi, the areas of the individual lung lobes dis- 
playing lesions were scored. Lung lobes revealed varying degrees 
of gross lesions, ranging from 0 to 75% (Fig. S2). Histological 
lesions in MERS-CoV-infected macaques were limited to the 
lungs. On day 3 postinfection, lesions were characterized as 
multifocal, mild (Fig. 54) to marked (Fig. 5B) interstitial pneu- 
monia. The pneumonia was characterized by thickening of al- 
veolar septae by edema fluid and fibrin, and small to moderate 
numbers of macrophages and fewer neutrophils. Alveoli contained 
moderate numbers of pulmonary macrophages and neutrophils; 
many of these alveolar macrophages formed multinucleate 
syncytia (Fig. 5D). One animal euthanized 3 dpi also displayed lesions 
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Fig. 2. Virus shedding in rhesus macaques inoculated with MERS-CoV. Nasal 
swabs (A) and BAL fluid (B) were collected 1, 3, and 6 dpi. Of note, CoV-2, 
CoV-4, and CoV-6 were euthanized 3 dpi and thus no samples are available 
for 6 dpi. RNA was extracted and viral load was determined as TCIDso 
equivalents by qRT-PCR. TCIDsq equivalents were extrapolated from stan- 
dard curves generated by adding dilutions of RNA extracted from a HCoV- 
EMC/2012 stock with known virus titer in parallel to each run. Animal num- 
bers are indicated on the x axis. Black bars indicate viral load at 1 dpi, white 
bars indicate viral load at 3 dpi, and gray bars indicate viral load at 6 dpi. 
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Fig. 3. Viral load in respiratory tissues of rhesus macaques inoculated with 
MERS-CoV. Rhesus macaques were euthanized on day 3 (black bars) and day 
6 (white bars) postinfection and tissue samples were collected. RNA was 
extracted and viral load was determined as TCIDsq equivalents by qRT-PCR. 
TCIDsq equivalents were extrapolated from standard curves generated by 
adding dilutions of RNA extracted from a HCoV-EMC/2012 stock with known 
virus titer in parallel to each run. Geometric mean viral loads were calcu- 
lated; error bars represent SD. R, right; L, left; LN, lymph node. 


similar to those described as bronchiolitis obliterans organizing 
pneumonia (Fig. 5C). This organizing pneumonia consisted of 
ageregations of fibrin (Fig. $3), macrophages, and sloughed pul- 
monary epithelium that fill and occlude small airways, an early 
lesion of bronchiolitis obliterans organizing pneumonia that 
has yet to develop the characteristic dense collagenous stroma, 
and marked infiltration of fibroblasts that would occur if the 
lesion was given time to mature. Perivascular infiltrates of in- 
flammatory cells were observed multifocally within and adjacent 
to affected areas of the lung. Animals necropsied at 6 dpi dem- 
onstrated moderate (Fig. 64) to marked (Fig. 6B) changes in 
the lungs with type II pneumocyte hyperplasia (Fig. 6C), 
abundant alveolar edema, and fibrin with formation of hyaline 
membranes (Fig. 6D). 


MERS-CoV Replicates in Alveolar Pneumocytes. At 3 dpi, in situ 
hybridization (ISH) demonstrated viral RNA multifocally in the 
lungs of all monkeys. The signal was found in spindle to polyg- 
onal cells along alveolar septae (Fig. 5). This location and the 
cellular morphology suggested that virus replication occurs 
mainly in the pulmonary epithelium. Viral RNA was not detec- 
ted in endothelial cells. Immunohistochemistry (IHC) yielded 
similar results for the presence of viral antigen (Fig. 5F). Addi- 
tionally, there were rare round mononuclear cells and stellate 
cells within the cortex of the mediastinal lymph nodes that were 
positive for both viral RNA by ISH and viral antigen by IHC. 
These cells were morphologically consistent with cells of the 
macrophage or dendritic cell lineage. Viral RNA and viral an- 
tigen could still be detected at 6 dpi (Fig. 6 F and F). 

To confirm that virus replication occurred mainly in alveolar 
pneumocytes, fluorescent IHC was performed using anticytokeratin 
and anticoronavirus antibodies. Viral antigen was predomi- 
nantly found in alveolar pneumocytes on 3 as well as 6 dpi (Fig. S4 
A and B). This distribution of viral antigen in type I and type II 
pneumocytes corresponded with the distribution of the receptor 
for MERS-CoV, DPP4 (Fig. S4 C and D) as indicated by IHC 
performed on the lungs of infected animals. 
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Analysis of Gene Expression in Lungs and PBMC. To assess genes that 
may play a significant role in disease pathology, we performed 
a microarray analysis on lung tissue and PBMC. A Welch’s ¢ test 
(P < 0.01, fold-change > 2) was used to compare profiles of 
significant differentially expressed genes (DEG) from the right 
lower lung lobe to observable gross lesions, relative to uninfected 
controls. The right lower lobe was selected, as qPCR demon- 
strated that viral genomes were present in those samples, 
allowing us to compare sites of histopathologic injury to the 
background transcriptional profiles of sites exposed to virus. On 
day 3, we observed a total of 173 DEG (Fig. 7A and Table 82), 
although the expression was highly variable between animals and 
sample sites. This variability is likely multifactorial, driven by 
differences between individual animals, such as levels of virus in 
the lung and at lesion sites, distribution of virus receptors, lung 
sizes, age, and distribution and quantity of susceptible cells. 
Using ingenuity pathway analysis (IPA), we determined that the 
majority of DEG were associated with antiviral immunity, in- 
flammation, and chemotaxis. Notably, genes such as IL-6, che- 
mokine C-X-C ligand 1 (CXCL1), and matrix metalloproteinase 
9 (MMP%9) associated with proinflammatory processes and re- 
cruitment of inflammatory cells were significantly up-regulated in 
the lesions compared with the right lower lung lobes. By day 6, 
we identified only 37 DEG distinguishing right lower lung lobes 
from lesions (Fig. 7B and Table S3), and these were not grouped 
in any heavily enriched functional categories. Additionally, by 
6 dpi the lung lesions showed a higher magnitude of gene ex- 
pression compared with the right lower lung lobe samples. This 
finding may reflect the resolving clinical status of the animals at 
day 6 postinfection, and indicate that substantial differences in 
gene expression occur at sites of virus-induced injury rather than 
throughout the lung. 

We sought to determine if MERS-CoV-induced gene ex- 
pression dynamics can be observed in PBMC and lung tissue and 
lesion samples from infected macaques, and change over time. 
Using singular value decomposition-coupled multidimensional 
scaling, we identified changes in gene expression occurring in 
parallel in both organs (Fig. S5). In PBMC a significant response 
to MERS-CoV infection is observed only on day 1 postinfection. 
In lungs a significant global response was found only 3 dpi (Fig. 
S5). By 6 dpi, the gene-expression dynamics in both organs 
return close to an uninfected or baseline state. Furthermore, 
functional analysis with IPA reveals that MERS-CoV infection in 
macaques leads to very rapid innate immune activation (14) 
through pattern-recognition receptors. Leukocyte activation is 
evident in blood and lungs, as well as evidence of a robust but 
self-limiting inflammatory response. Importantly, innate immune 


Fig. 4. Virus replication in the lungs of MERS-CoV inoculated rhesus mac- 
aques. Transmission electron microscopy was performed on lung lesion 
samples collected from MERS-CoV-inoculated macaques at 6 dpi. Virus 
particles (arrows) are visible in type | pneumocytes (Left and Right). (Scale 
bar, Right, 250 nm; Inset, 50 nm). B, basement membrane; E, endothelial cell; 
P, pneumocytes; R, red blood cell. 
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Fig. 5. Histopathological changes at 3 dpi in rhesus macaques inoculated 
with MERS-CoV. Rhesus macaques were euthanized on day 3 postinfection 
and lung tissue was collected and stained with H&E. (A) Mild acute in- 
terstitial pneumonia demonstrated by mild thickening of alveolar septae by 
small numbers of neutrophils, macrophages, fibrin, and edema. Small 
numbers of inflammatory cells and fibrin also extend into alveolar spaces. (B) 
Marked interstitial pneumonia with bronchiolitis obliterans (arrow). Airways 
are filled with organizing fibrin, necrotic debris, collagen and fibroblasts. 
There are numerous multinucleate syncytial cells scattered throughout the 
alveoli (arrowhead). (C) Bronchiolitis obliterans organizing pneumonia. (D) 
Alveolar capillary (asterisk) bounded by degenerate, vacuolated type | pul- 
monary epithelium. The alveolar space contains inflammatory cells and 
a multinucleate syncytial cell (arrowhead). (E) ISH demonstrates viral RNA 
predominantly within alveolar pneumocytes. (F) Anticoronaviral IHC reveals 
viral antigen within alveolar pneumocytes. (Magnification: A and B, 200x; C, 
E, and F, 400x; D, 1,000x.) 


activation occurs very rapidly after infection, and is subsequently 
also rapidly controlled, indicating quick resolution. 


Cytokine and Chemokine Profiles in MERS-CoV-Infected Animals. 
Plasma levels of 23 cytokines and chemokines were determined 
in serum samples obtained at 0, 1, 3, and 6 dpi. Transient dif- 
ferences were only observed for IL-1RA, IL-2, IL-6, IL-8, IL-12, 
IL-13, IL-15, IFN-y, and monocyte chemotactic protein-1 (MCP- 
1). AT 1 dpi serum levels of IL-1RA, IL-2, IL-13, IL-15, INF-y, 
and MCP-1 were significantly higher than at 0 dpi and all levels 
returned to baseline by 6 dpi (Fig. S6). Although the serum level 
of IL-6 was higher and serum levels of IL-8 and IL-12 were lower 
at 1 dpi, these differences were not statistically significant. 


Discussion 


The emergence of MERS-CoV in 2012 was the second in- 
troduction in 10 y of a novel coronavirus, causing severe re- 
spiratory disease, into the human population. Both the severity 
of the respiratory disease and the genetic similarity to betacor- 
onaviruses isolated from bats are reminiscent of the emergence 
of SARS-CoV in 2003, which resulted in 8,422 human cases with 
916 deaths in over 30 countries (15). To understand the disease 
caused by MERS-CoV in humans, we developed a nonhuman 
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primate model, which recapitulates the respiratory disease ob- 
served in mild or moderate human cases but not the occasionally 
occurring kidney disease. 

In this rhesus macaque model, virus shedding as indicated by 
qRT-PCR occurred predominantly via the nose and, to a limited 
extent, the throat. In nasal swabs and BAL, viral loads were 
highest day 1 postinfection and decreased over time. However, at 
6 dpi, two of three animals were still shedding virus from the re- 
spiratory tract. Although MERS-CoV was detected in the upper 
respiratory tract and the lymphoid tissue draining the lungs, 
replication of MERS-CoV was most prominent in the lower 
respiratory tract. MERS-CoV replicated predominantly in type 
I and II pneumocytes in the alveoli. These two cell types form 
the main component of the architecture of the alveolar space 
around the terminal bronchioles. This predominant replication 
of MERS-CoV in alveoli may explain the limited amount of 
virus shedding observed in our rhesus macaque model. In ad- 
dition, the fact that human-to-human transmission so far seems 
limited to a few family clusters in Saudi Arabia (6), the United 
Kingdom (7), and Tunisia and nosocomial transmission in Jordan, 
Saudi Arabia (6), and France (9) might be explained by the pro- 
pensity of the novel coronavirus to replicate deep down in the 
lower respiratory tract. Recently, DPP4 was identified as the re- 
ceptor for MERS-CoV (12). DPP4 is an evolutionarily conserved 
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Fig. 6. Histopathological changes at 6 dpi in rhesus macaques inoculated 
with MERS-CoV. Rhesus macaques were euthanized at 6 dpi and lung tissue 
was collected and stained with H&E. (A) Mild acute interstitial pneumonia 
demonstrated by mild thickening of alveolar septae by small numbers of 
neutrophils, macrophages, fibrin, and edema. Small numbers of inflammatory 
cells and fibrin also extend into alveolar spaces. (B) Marked interstitial pneu- 
monia with extensive type Il pneumocyte hyperplasia (arrow). Alveolar spaces 
are flooded with edema fluid, fibrin, and inflammatory cells. (C) Type II 
pneumocyte hyperplasia (arrow) with alveolar flooding edema and fibrin. 
(D) Hyaline membranes line alveolar space (arrowhead). (E) ISH demonstrates 
viral RNA predominantly within alveolar pneumocytes. (F) Anticoronaviral IHC 
reveals viral antigen within alveolar pneumocytes. (Magnification: A and B, 
200x; C-F, 400x.) 
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Fig. 7. Microarray analysis of gene expression in lungs of rhesus macaques 
inoculated with MERS-CoV. The heatmap shows logj9 ratios to uninfected 
controls of differentially expressed genes in the right lung lower lobe and 
lung lesion samples as determined by Welch’s t test (P < 0.01, fold-change > 
2), and grouped by hierarchical clustering. A representation is shown of 173 
differentially expressed genes from individual animals and sample sites at 3 
dpi (A) and of 37 differentially expressed genes from individual animals and 
sample sites at 6 dpi (B). Functional analysis was performed on subgroups of 
the day 3 differentially expressed genes, and enriched functional categories 
identified by IPA are shown in the blue and red vertical bars to the right of 
the heatmap in A. Differentially expressed genes are specified in Table S2 
(3 dpi) and Table S3 (6 dpi). 


type II transmembrane glycoprotein expressed on endothelial and 
epithelial cells, in lungs and kidney. Currently, the exact distribu- 
tion of DPP4 throughout the human respiratory tract is unknown, 
although DPP4 expression has been observed in vitro in human 
bronchiolar epithelial cells (12). In mouse lungs, DPP4 is present 
on type I and type II pneumocytes in the alveoli, with expression 
on type I pneumocytes more abundant than on type II pneumo- 
cytes (16). The distribution of DPP4 is in agreement with our 
observation of MERS-CoV virus replication in the alveoli and 
suggests a role for DPP4 in tissue tropism of the novel coro- 
navirus. A more thorough analysis of DPP4 distribution in the 
human airway is needed to examine the apparent alveolar re- 
striction of virus replication and whether efficient human-to- 
human transmission might be restricted by the receptor distribution, 
similar to what has been described for highly pathogenic avian 
influenza A/HSN1 viruses (17, 18). No virus replication was ob- 
served in kidney samples from MERS-CoV-infected rhesus mac- 
aques. This finding could indicate that the renal failure reported in 
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some of the human cases might be a result of hypoxia rather than 
virus replication in the kidney, or of other underlying, preexisting 
medical conditions. 

Currently, few clinical and no pathological data are available 
on the human cases of MERS-CoV. In addition, the data that 
are available might be biased toward severe cases in need of 
hospitalization, as there is no information on the prevalence of 
MERS-CoV in the general population. Recent reports have 
revealed that the majority of the severe human MERS-CoV 
cases have underlying medical conditions, potentially negatively 
affecting outcome of disease in these patients (2, 8, 9, 19). As 
such, a direct comparison of the pathogenesis of the novel 
coronavirus observed in the rhesus macaque model to that in 
humans should be approached cautiously. However, the tran- 
sient nature of disease observed in rhesus macaques suggests that 
this model does not exactly recapitulate the disease observed in 
the fatal cases of MERS-CoV, likely because of the health status 
of these animals compared with severe human MERS cases. 
Rather, the MERS-CoV-associated disease observed in rhesus 
macaques, and the variation observed between individual ani- 
mals, may be more similar to the recent, less-severe cases of 
MERS-CoV identified in the United Kingdom and Saudi Arabia. 
Interestingly, the disease displayed in the rhesus macaques is 
reminiscent of SARS-CoV in nonhuman primates, because a wide 
disease spectrum was reported in the different SARS-CoV studies 
in nonhuman primates (20-27). 

In the case of SARS-CoV, the ability to cause severe disease in 
a macaque model was associated with an exacerbated innate host 
response in aged animals (10-19 y old) (28). The majority of 
human cases of MERS-CoV have been middle aged and aged 
male individuals (39-73 y old) (6). In the rhesus macaques in- 
oculated with MERS-CoV, no difference in severity of disease 
could be attributed to age (6-12 y old, young adult) or sex, but 
low animal numbers used here do not allow a firm conclusion 
and we cannot exclude that severe disease could be age de- 
pendent in the MERS-CoV model described herein. 

Proinflammatory cytokines are likely to play an important role 
in the observed acute lung injury observed in the fatal human 
cases. The up-regulation of the expression levels of proin- 
flammatory cytokines and chemokines, such as IL-6, CXCL1, 
and MMP%9 in the rhesus macaques and the associated cytokine 
and chemokine levels in serum samples may have resulted in 
chemotaxis and activation of neutrophils. Indeed, increased 
numbers of neutrophils were observed in the blood of all animals 
on 1 dpi, in accordance with what has been reported for the first 
human case (2). 

The transient, self-limiting nature of MERS-CoV infection in 
rhesus macaques was displayed at various levels within the 
model. Transient patterns were observed for virus shedding, viral 
loads in tissues, neutrophilia, type Il pneumocyte hyperplasia, 
regulation of host gene expression, and production of cytokines. 
The controlled IFN-mediated innate immune response was very 
rapid and resolved within the first 3 d. Correlation of gene- 
expression dynamics between PBMC and lung-tissue samples also 
suggests that circulating profiles may have utility as a means of 
assessing host responses to infection without requiring an invasive 
tissue biopsy. 

Taken together, our study provides data on the extent and 
duration of respiratory disease caused by the novel coronavirus 
in rhesus macaques. MERS-CoV primarily affects the epithelium 
of the lower respiratory tract in a transient manner. This tropism 
of MERS-CoV for the lower respiratory tract explains both the 
severity of the disease observed in humans and the limited 
number of cases of human-to-human transmission to date. 


Materials and Methods 


Ethics Statement. All animal experiments were approved by the Institutional 
Animal Care and Use Committee of the Rocky Mountain Laboratories, and 
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performed following the guidelines of the Association for Assessment and 
Accreditation of Laboratory Animal Care, International (AAALAC) by certified 
staff in an AAALAC-approved facility. 

The work with infectious MERS-CoV was approved under biosafety level- 
3 conditions by the Institutional Biosafety Committee. Sample inactivation 
was performed according to standard operating procedures approved by 
the Institutional Biosafety Committee for removal of specimens from high 
containment. 


Virus and Cells. MERS-CoV isolate HCoV-EMC/2012 was kindly provided by the 
Viroscience Laboratory, Erasmus Medical Center, Rotterdam, The Nether- 
lands, and propagated in VeroE6 cells in DMEM (Sigma) supplemented with 
2% (vol/vol) FCS (Logan), 1 mM L-glutamine (Lonza), 50 U/mL penicillin, and 
50 pg/mL streptomycin (Gibco). 


Animal Study and Sample Collection. Six rhesus macaques (four males, two 
females; age 6-10 y, weight 5-10 kg) were inoculated with 7 x 10° TCIDso of 
MERS-CoV isolate HCoV-EMC/2012, as described previously (13). For detailed 
information, see S/ Materials and Methods. 


Histopathology, ISH, and IHC. Histopathology, ISH, and IHC were performed on 
macaque tissues. For detailed information, see S/ Materials and Methods. 


Transmission Electron Microscopy. Transmission Electron Microscopy was per- 
formed essentially as described previously (29). For detailed information, see 
SI Materials and Methods. 


RNA Extraction. RNA was extracted from swabs, BAL, and whole blood 
samples using the QiaAmp Viral RNA kit (Qiagen) according to the manu- 
facturer’s instructions. Tissues (30 mg) were homogenized in RLT buffer and 
RNA was extracted using the RNeasy kit (Qiagen). 
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Quantitative PCR. For detection of viral RNA, we used the one-step real-time 
RT-PCR upE assay (30); for the detection of mRNA, we used an ORF1- and 
ORF5-specific SYBR green assay. See S/ Materials and Methods for detailed 
information. 


Microarray Data and Functional Analysis. RNA extracted from lung tissue 
samples and PBMC was run on Agilent Rhesus macaque 4 x 44 K microarrays 
(Agilent Technologies) using the manufacturer's one-color analysis protocol. 
Raw array data were uploaded to Genedata Analyst 7.6 (Genedata) and 
further analyzed. Raw microarray data have been deposited in National 
Center for Biotechnology Information’s Gene Expression Omnibus (GEO) 
database under GEO accession no. GSE44542, and are also available to the 
public at http://viromics.washington.edu. For detailed information, see S/ 
Materials and Methods. 


Serum Cytokine and Chemokine Analysis. Concentrations of cytokines and 
chemokines in serum were measured using the Non-Human Primate Cytokine 
MILLIPLEX map 23-plex kit (Millipore) according to the manufacturer’s 
instructions. See S/ Materials and Methods for more information. 
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